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ABSTRACT 

As the number of known exoplanets continues to grow, the question as to whether such 
bodies harbour satellite systems has become one of increasing interest. In this paper, 
we explore the transit timing effects that should be detectable due to an exomoon 
and predict a new observable. We first consider transit time variation (TTV), where 
we update the model to include the effects of orbital eccentricity. We draw two key 
conclusions: 

1) In order to maintain Hill stability, the orbital frequency of the exomoon will always 
be higher than the sampling frequency. Therefore, the period of the exomoon cannot 
be reliably determined from TTV, only a set of harmonic frequencies. 

2) The TTV amplitude is cx Mgas where Ms is the exomoon mass and as is the semi- 
major axis of the moon's orbit. Therefore, Ms & as cannot be separately determined. 

We go on to predict a new observable due to exomoons - transit duration variation 
(TDV). We derive the TDV amplitude and conclude that its amplitude is not only 
detectable, but the TDV signal will provide two robust advantages: 

1) The TDV amplitude is oc M s a~ s 1/2 and therefore the ratio of TDV to TTV allows 
for Ms and as to be separately determined. 

2) TDV has a ir/2 phase difference to the TTV signal, making it an excellent comple- 
mentary technique. 

Key words: techniques: photometric — planets and satellites: general — planetary 
systems — occultations — methods: analytical 



1 INTRODUCTION 

Over 300 exoplanets have been discovered to date with 
detection rates escalating (see http: / / exo planet.euj by J. 
Schneider). The detections have been so far biased towards 
large bodies and the smallest transit ing planet to date is 
still Neptune-sized, for Gliese 436b dGillon et all l |2007ft h 
Current instruments cannot yet detect transiting Earths but 
transit time variation (TTV) could offer a way to bring sen- 
sitivit y do wn to sub Earth-m ass level (|Holman fc Murray! 
l |2005j ) and lAgol etafl (|2005h ). Given the large number of 
moons in our own Solar System, it is reasonable to postu- 
late that satellites are common around exoplanets. 

Another profound motivation for looking for exomoons 
is that they are likely to be terrestrial in nature, based upon 
our own Solar System, and hence one would propose that 
exomoons could be more habitable environments than the 
host of extrasolar giant planets (EGPs) so far discovered. 

Photometric detection of a moon is likely to be excep- 
tionally challenging. An exomoon is likely to be sub-Earth 
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sized, based upon our Solar System, and so even a special- 
ized transit space-base d telescope like COROT will st ruggle 
to spot the signature |Sartoretti fc Schneider] |l999)). An 
additional problem lies in the fact that much of the time 
the exomoon will not appear to be its orbital distance from 
the planet, but some fraction of it, depending on the orbital 
phase of the exomoon during transit. The moon can effec- 
tively hide behind the planet or in front of it. This makes 
disentangling the photo metric signature excepti o nally diffi- 
cult and was outlined by Sartoretti fc Schneider! (|l999h and 
ICabrera fc Schneider! \200ij ). 

Pre viously, several differ en t au thors ( Sza bo et al 

(2006), ISartoretti fc Schneider! ||l999h and I imon et al 



(2007)) have noted that transit timing variations (TTV) 
could be used to indirectly detect the presence of such exo- 
moons. In the hrst half of this paper, we update the model 
for the TTV effect by including orbital eccentricity. We com- 
pare our formulation to the previously proposed mathemat- 
ical treatment and demonstrate our model reduces to the 
original analytic equations for zero eccentricity. 

However, the crucial problem with TTV is that the 
amplitude of the signal is proportional t o both the exo- 
moon mass, Ms, and orbital separation, q.<?. lFord fc Holmanl 
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(120071 1 referred to this obstacle as the 'inverse problem', in 
regard to transit timing effects. In §3.1 we show that the 
exomoon's period cannot be reliably determined from TTV 
and hence as remains an unknown. Ergo, one cannot es- 
tablish the mass of the exomoon without assuming a value 
for the orbital separation. It is therefore clear that a strong 
desideratum is a secondary method which can complement 
TTV and remove this degeneracy and this constitutes the 
focus of the second part of our paper. 

By considering the transit duration, we predict a new 
observable timing effect due to an exomoon, which we label 
as transit duration variation (TDV) . We find that the ampli- 
tude of this timing signal is of the same order of magnitude 
to the TTV signal and indeed often larger. The effect is also 

\ j 2 

predicted to be oc Msa s . Hence, the ratio of TDV to 
TTV allows for the mass of the exomoon to be found with- 
out assuming an orbital distance. In addition, TDV is 7r/2 
out-of-phase with TTV, making it an ideal complementary 
method for exomoon detection. 



2 TTV AMPLITUDE DUE TO AN EXOMOON 
2.1 Outline of the model 

In the first half of this paper, we aim to update the model for 
the transit timing variation (TTV) signal due a transiting 
exoplanet with a single satellite of mass Ms and non-zero 
orbital eccentricity. In this work, we consider the variation 
of the mid-tran sit point of th e plan etary transit, Tmid - This 
is in contrast to lSimon et al.l (|2007n who consider the planet 
and moon combined transit. Throughout our discussion, we 
also make the assumption that planet-moon orbital plane is 
co-aligned with the planet-star orbital plane at i = 90°. 

In our case, the planet orbits the barycentre of the 
planet-moon system with a semi-major axis of aw where 
W denotes wobble. The fact that aw > means that the 
time between the planet being at the mid-transit point and 
the barycentre being at the mid-transit point is, in general, 
non-zero, and this is the origin of the TTV effect. 

Consider the projected distance between the planet and 
the planet-moon barycentre to be given by x' 2 , as illustrated 
in figure 1. The TTV effect will be given by x' 2 divided by the 
^-direction component of the barycentre's orbital velocity 
around the star, given by vb±- Since x' 2 is a function of the 
planet's true anomaly around the planet-moon barycentre, 
fw, so too is the TTV effect: 



TTV(/ V 



x' 2 (fw 

VB± 



(1) 



In the case of a circular orbit, we expect the TTV sig- 
nal to have a sinusoidal nature, and infact the peak-to-peak 
amplitude of the wave may be unambiguously defined as 
2aw /vb±- In the case of eccentric orbits, we expect non- 
sinusoidal waveforms, due to Kepler's Equation. Peak-to- 
peak amplitudes can become ambiguous in such cases, and 
so we choose to use the root-mean-square (rms) amplitude 
definition, which is valid for all waveforms. For a simple cir- 
cular orbit, the rms amplitude of the TTV signal, Sttv, will 
be given by: 




Figure 1. Cartoon of the star-planet-moon system. In this 
schematic, the observer lies at y 2 = +oo and the exomoon is not 
shown, but the wobble of the planet due its presence is represented 
by the smaller ellipse. The origin of the x' 2 -y 2 scale denotes the 
position of the planet-moon barycentre, which itself orbits the star 
in the x 2 -y 2 scale. 
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Where the distance aw is small enough that we may 
assume vb± is a constant over the time-scale of the TTV 
effect. 

Let us now extend the analysis to include the effect of 
exomoon orbital eccentricity, eg and position of pericentre, 
zus, as well as planetary eccentricity ep and position of peri- 
centre, zup. In appendix A, we derive the general equation 
for the rms amplitude of the TTV effect due to an exomoon 
to be given by: 



(rjes, tgs) 
V2 y/G(M» + Mprv) T(e P ,zu P ) 



(3) 
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T = cos arctan 
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Where ap is the semi-major axis of the exoplanet's orbit 
around the central star, as is the semi-major axis of the 
exomoon's orbit around the exoplanet, G is the gravitational 
constant, M* is the mass of the host star and Mprv = Ms + 
Mp i.e. the combined mass of the exomoon and exoplanet 
respectively. 



Mppy is the mass of the planet as measured by radial velocity 
surveys, and thus includes the mass of satellites 
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In the case of es — * we have (j- -t 1 and simi- 
larly for ep — > we have T — > 1. Thus we may compare 
our equation to that derived by previous authors, who as- 
sume d circular orbits, and expect to produce the same re- 
sult. ISartoretti fc Schneider! f|l999l ') predicted the peak-to- 
peak amplitude of an exomoon around an exoplanet with 
orbital period Pp to be given by equation (6). 



At 3 



2a s M s Mp 1 x P P (2ivap) 



(6) 



Using Kepler's Third Law, it is trivial to show that 
this is equivalent to our expression, for circular orbits and 
Mp 3> Ms, except for a factor of 2V2, which comes from 
the fact we use an rms definition of amplitude, rather than 
a peak-to-peak definition. Thus we have confirmed that 
our general expression for the TTV effect due to an exo- 
moon is equivalent to the circ ular equations first derived by 
ISartoretti fe Schneider! ([1999). For completion, we note that 
the TTV amplitude may be written purely in terms of the 
masses of the three bodies and the period of the exoplanet 
(see appendix A, equation (A28) & (A29)). 

The effect of eccentricity is discussed in more detail in 
§3.3 and typical waveforms produced are shown later in fig- 
ure 3, §4.2. In §4.2, Table 1 gives a list of the expected TTV 
rms amplitudes due to a liWfg exomoon around some of the 
most promising exoplanet candidates. 



2.2 Permitted range for as 

The orbital radius of any satellite around a planet must lie 
somewhere between the Hill radius, d ma x, and the Roche 
limit, d m in, to maintain stabilitjjf). We express as by assum- 
ing that it is equal to some fraction, \, of the Hill radius, 

dmax - 



as = X ■ 



dm 
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= Ri 
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Where ap is the semi-major axis of the planet around 
the host star, Rp is the planetary radius, \ is some real 
number between and 1, and p denotes density. 

X may be f urther constrained by noting that 
Barnes fe O'Brien! <|2002h estimated x i$ 1/3 and 
Domingos et al.l ( 20061 ) estimated x i; 1/2- This is because 
the Hill sphere is just an approximation, and in reality other 
effects, like radiation pressure or the Yarkovsky effect, can 
perturb a body outside of the sphere. We choose to use the 
conservative choice of x ^ 1/3 an d combining this limit with 
the Roche limit, and rewriting in terms of planetary period, 
Pp for M* Mp, we can estimate: 



models of the internal structure for such bodies. If we use 
Rs oc Ms 27 , then we can estimate ps ~ O.lMg 19 . 
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3 IMPLICATIONS 

3.1 The high frequency nature of an exomoon's 
TTV 

If we make the approximation that M* ^> Mp^y and 
Mprv 3> Ms, and employ Kepler's Third Law, we may 
provide an estimate for the ratio of the exomoon to planet 
orbital period, which can be shown (see appendix B) to be: 



Ps_ 
Pp 



(12) 



Since x ^ 1 we will always be in the regime where Ps < 
Pp. Taking \ ~ 1/3 gives a rough estimate of Ps/ Pp ~ 1/9 
and so the frequency of the signal is vs ~ 9Pp X . How- 
ever, the Nyquist frequency will be given by one half of the 
sampling rate (O.bPp 1 ), which represents the maximum fre- 
quency we can resolve without aliasing. 

The usual technique for detecting signals within data is 
to employ a periodogram, but this method will suffer from 
aliasing in the exomoon case. Therefore, we can only derive 
a set of harmonic frequencies which the exomoon's orbit 
could exhibit. We therefore conclude Ps cannot be reliably 
determined from the TTV signal. 



3.2 The limitation of TTV in determining M s 

As originally pointed by ISartoretti fe Schneider! (1 19991 ). the 
major problem with TTV is that one cannot determine the 
mass of the exomoon without making an assumption on the 
distance at which the moon orbits the planet. Using equation 
(3), it is possible to write the TTV amplitude as a function 
of purely the exomoon properties: 



5ttv oc Ms as 



(13) 



Therefore we can effectively only determine the moment 
of the exomoon. If one knows the period of the exomoon, 
then it is trivial to derive as using Kepler's Third Law, but 
as seen in §3.1, Pp cannot be reliably ascertained. 

This crucial limitation of exomoon TTV makes mass 
estimation unfeasible and the best we can ever do is merely 
provide evidence for an exomoon within a large mass range. 
This severe limitation will be resolved later in this paper. 



IValencia et al. (2006) offer a way of rewriting ps in 
terms of just the mass of the exomoon by adopting terrestial 

2 Note, we use the rigid body Roche limit for simplicity 



3.3 The effect of eccentricity 

The effects of orbital eccentricity on the TTV amplitude 
are all absorbed into the two parameters, Cr(es,^s) and 
T(ep,u?p). If we increase es from zero to unity, regardless 
of what value zus takes, £t decreases below 1 and hence the 
TTV amplitude will always decrease. This implies exomoon 
detections are biased towards satellites on circular orbits. 
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Figure 2. Behaviour o/T — 1 as a function of ep andvjp, which 
is effectively the factor by which the TTV is enhanced by for 
ep > 0. Positions of pericentre near 270° produce consistently 
enhanced TTV amplitudes. The flat line corresponds to ep = 
and each progressively larger amplitude wave is for ep = 0.3, 0.6 
and 0.9 respectively. 



The situation is more complicated for ep, where a non- 
zero ep makes the TTV amplitude significantly increased or 
decreased depending on zap. In figure 2, we plot T _1 as a 
function of cap for several different eccentricities and find 
that for ep > 0, the most favourable position of pericentre 
is wp ~ Sn/2. However, we point out a recent study by 
iKane fc von Braunl (|2008T l which predicts such exoplanets 
to possess a low transit probability. 

For the eccentric transiting exoplanets GJ436b, XO- 
3b, HAT-P-2b and HD17156b, T" 1 takes values 1.01, 1.03, 
0.94 and 0.47 respectively. However, the stability of satellites 
around eccentric exoplanets remains unclear. 



4 TRANSIT DURATION VARIATION (TDV) 

4.1 The TDV due to an exomoon 

Transit duration variation (TDV) is the periodic change in 
the duration of a transit (r) over many measurements, where 
we define r as the time between the 1 st and 4 th contact 
points. It has previo usly been discussed as a possible test of 
general relativity bv lPal fc Kocsisl l|2008l ). In this discussion, 
we consider the TDV due to an exomoon and conclude that 
it should produce a detectable signal. 

We first consider that the duration of a transit is in- 
versely proportional to the projected velocity of the planet 
across the star, vp±, and make the following assumptions: 

■ The projected velocity of the planet-moon barycentre 
during the transit, relative to the star, does not vary signif- 
icantly over the time-scale of the transit duration. 

■ The projected velocity of the planet during the transit, 
relative to planet-moon barycentre, does not vary signifi- 
cantly over the time-scale of the transit duration. 

■ The orbital inclinations of the planet (i = 90°) and 
exomoon do not vary from orbit to orbit. 

■ We do not consider additional perturbing bodies in the 
system. 



t oc l/vpx (14) 

In this case, any transit duration variation must be 
solely due to the variation of the velocity. For a single com- 
panion exomoon, the velocity of the planet will have two 
components: 

dpi = v B± + v W x (15) 

where vp. is the velocity of the planet-moon barycentre 
around the host star and vw is the velocity of the planet 
around the planet-moon barycentre, i.e. the wobble of the 
planet due its companion satellite. 

It is clear that vw± will be significantly different for 
each transit unless Pp / Ps is some low-order integer. Some- 
times vw± will be additive to the barycentre's velocity and 
sometimes subtractive resulting in shorter and longer transit 
durations respectively; an effect we label TDV. 

The factor by which r will vary must be equal to the 
ratio of the velocities R — vw±/vb± and from this starting 
point the TDV amplitude may be shown (see appendix C2) 
to be equal to: 

s _ fap / Ms _L Co( e s.^s) 

TDV -\I as ' y Mprv(M prv + M») ' y/2 ' T(ep,^p) 

(16) 

where 

j. , , /l + e% — e% cos(2ros) 
Ci>(es,^s) = W s (17) 

In an analogous way to the TTV amplitude, this may be 
re-written in terms of simply the masses in the system, and 
this equation may be found in the appendix C2, equations 
(C28) and (C29). We also point out that this effect has the 
following proportionality: 

Stdv oc M s ■ a~ 1/2 (18) 

As a final note, we point out that in §3.3 we discussed 
how increasing eg tends to decrease the TTV amplitude. For 
the TDV signal, the opposite is true, increasing es tends to 
increase the TDV signal. 

4.2 TTV & TDV as complementary methods 

From equations (13) and (18), it is clear that the ratio of 
TDV to TTV should be able to eliminate Ms and we can di- 
rectly measure the orbital seperation of the exomoon as and 
hence Ms- The introduction of TDV allows for the precise 
measurement of Ms without any assumption on the orbital 
separation. Using Kepler's Third Law, it is then possible to 
derive the exomoon's orbital period, which was shown in 
§3.1 to be unattainable from TTV alone. In appendix D, it 
is shown that if we assume es — (but ep can take any 
value between zero and unity), then r\ = Stdv /Sttv is ap- 
proximately given by: 

OTDV 

71 ~ 5^r7 ~ Tp~ ' ^372 ( ' 
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Figure 3. The TTV and TDV waveforms due to a hypothetical 
exomoon around GJ436b. We show a) the effects of exo- 
moon eccentricity: the solid line gives eg = 0, the dashed line 
e s = 0.3 and the dotted line e s = 0.6 b) TTV leads TDV by a 
■k/2 phase difference c) the high frequency nature of both waves 
(as discussed in §3.1) is apparent d) increasing eg tends to de- 
crease the TTV amplitude and increase the TDV amplitude (as 
discussed in §3.3 & §4.B) 



Therefore one may determine % an d hence as (and Ms) 
if one assumes the exomoon has zero ecce ntricity. This as- 
sumpt ion is supported based on a study bv lDomingos et all 
(2006) and the pattern observed amongst large (> 500km) 
moons in our own Solar System. 

Another major advantage of TDV is that the signal 
should lag TTV by a n/2 phase difference, originating from 
the fact TTV is a spatial effect whereas TDV is a velocity 
effect. Unfortunately, the direction and value of the phase 
shift remains unchanged between prograde and retrograde 
satellites. 

Combining TTV and TDV should allow for a much 
more significant confirmation of a potential exomoon than 
just using TTV alone. The phase difference can be seen in 
the example waveforms in figure 3, where we plot the TTV 
and TDV waveforms due to a hypothetical exomoon of IMq 
around a GJ436b for es = 0, 0.3 and 0.6. 

In Table 1, we predict the TTV and TDV rms ampli- 
tudes for a variety of known transiting planets. We include 
the effects of the known planetary eccentricity, and consider 
a lMgj exomoon with x = 1/3 and es = 0. 

In most cases, the TDV amplitude is actually greater 
than the TTV amplitude. We do however point out that the 
error in the transit durati on is typically 2-3 t imes that of 
the mid-transit time (e.g. lAlonso et al, r <|2008f >). Neverthe- 
less, the signal should be detectable, especially with future 
telescopes like Keplei0. It is interesting to note that timing 
errors of just under 1 second could detect a lM^ around 
a 11.79Mj exoplanet in the case of XO-3b, a mass ratio of 
~ 3 x 10~ 4 . 



4.3 A hypothetical habitable exomoon 

Consider a hypothetical exoplanet which is identical to 
GJ436b except it has a period of Pp = 35.7 days and is on 
a circular orbit, putting it into the habitable zone of GJ436 
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$TDV 1 


GJ436b 


14.12 


13.68 


CoRoT-Exo-4b 


7.58 


9.15 


OGLE-TR-lllb 


4.63 


7.32 


HAT-P-lb 


4.58 


6.82 


HD149026b 


3.61 


9.76 


Lupus-TR-3b 


3.28 


5.19 


WASP-7b 


3.26 


5.88 


TrES-lb 


3.04 


5.95 


HD17156b 


3.07 


1.06 


HD209458b 


2.97 


5.95 


XO-5b 


2.65 


4.69 


HAT-P-4b 


2.54 


8.34 


HD189733b 


1.52 


2.96 


XO-3b 


0.41 


0.87 



see http://kepler.nasa.gov/sci/ 



Table 1. Predicted TTV & TDV rms amplitudes due to a lAf® 
exomoon, for a selection of the best candidate transiting planets. 
System parameters are taken from various references, which are 
shown. 

Reference 



Alonso et al. (2008) 
Aigrain et al. (20081 

Di'az et al. (2008) 
Johnson et al. (20081 

Winn et al. (20071 
Weldrake et al. (2008) 
Hellier et al. (20081 
Winn et al. (2007a) 
Barbieri et al. (20071 
Kippine (2008) 
Burke et al. (20081 
Kovacs et al. (20071 
Winn et al. (2007b1 
Winn et al. (20081 



with T cq ~ 300K. Although the Neptune-like planet itself 
would not be an ideal place to search for life, an Earth mass 
exomoon would be. Suppose there is a 1M^ exomoon with 
this planet on a circular orbit and \ = 0-25 and hence a pe- 
riod of Ps — 2.5 days. Could we detect such a provocative 
target using timing effects alone? 

The predicted rms TTV amplitude woul d be 138s and 
the TD V amplitude would be 60s. For GJ436b. lAionso et all 
(2008) used the 1.52m Telescopio Carlos Sanchez telescope 
and achieved a timing accuracy of ~ 13 seconds and a dura- 
tion error of ~ 50 seconds. This suggests that the detection 
of the exomoon should be presently possible through TTV 
from the ground, and feasible with TDV in the near future. 
This illustrates that even ground-based instruments could 
detect an Earth-like body in the habitable zone using tim- 
ing effects. 



5 DISCUSSION & CONCLUSIONS 

We have presented an updated model for the TTV signal due 
to an exomoon to include the effects of orbital eccentricity 
in both the exoplanet and the exomoon. From the updated 
TTV model, we draw the following conclusions: 

(a) TTV is degenerated in that it can only determine 
Ms x as, where Ms is the exomoon mass and as is the 
exomoon's orbital radius. 

(b) The TTV due to an exomoon can be significantly en- 
hanced for exoplanets of ep > and zup ~ 270°. However, 
it remains unclear how dynamically stable such exomoons 
would be. 

(c) The TTV frequency will be always be greater than 
the sampling frequency of once every transit, implying we 
can only determine a set of possible harmonic frequencies 
for the exomoon's period, Ps- 

An exomoon is predicted to have another detectable 
timing effect on a transit in the form of transit duration 
variation (TDV). We have derived an equation for predicting 
the TDV amplitude and drawn the following conclusions: 
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(a) The ratio of TDV to TTV allows for the separate de- 
termination of Ms and as; thus solving the 'inverse problem' 
for exomoons. 

(b) The TDV signal is of a similar order of magnitude to 
the TTV signal. 

(c) the TDV signal is also enhanced for ep > and zap ~ 
270°. 

(d) TDV lags TTV by a 90° phase difference, making it an 
excellent complementary technique for exomoon detection. 

We also find that current ground-based telescopes could 
detect a lMgj exomoon in the habitable zone around a 
Neptune-like exoplanet. The author would therefore encour- 
age observers to produce not only their mid-transit times, 
but also transit durations for each transit, rather than com- 
posite lightcurve durations. This will allow constraints to be 
placed on the presence of exomoons around such planets. 
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APPENDIX A: TTV RMS AMPLITUDE 

The root-mean-square amplitude of a waveform with dis- 
placement TTV (fw) as a function of some variable fw may 
be written as: 



&TTV = \hr [TTV(/w)] 2 d/w 



(Al) 



where fw is true anomaly of the planet during its orbit 
around the planet-moon barycentre and is the only variable 
we integrate over and thus we are assuming: 

• ep & zap do not change over the measurement time- 
scale of the TTV effect. 

• There are no other secular changes in the planet's orbit. 

• es does not change over the measurement time-scale of 
the TTV effect. 

• zus takes the same value when sampled once every plan- 
etary orbital period, but fs and fw do not. 

In our case, the function TTV(fw) is given by the pro- 
jected distance between the planet and the planet-moon 
barycentre divided by the projected velocity of the barycen- 
tre, as a function of the planet's true anomaly around the 
planet-moon barycentre, fw- 



TTV(/ M 



x' 2 (fw,es,zus) 
vs.±(ep,ujp) 



(A2) 



Where x' 2 (fw, es,vus) is the projected displacement of 
the planet away from the planet-moon barycentre (see figure 
1), as a function of fw- 

To find x'ytfw . eg , zus), we take a similar approach to 
that of iKippina l|2008l ). where the author considered an ini- 
tial frame S' and then made a series of transformations. If 
we start out with the same setup as the cited author, we 
have an ellipse located with the centre at the origin of an 
x'-y' plot and the barycentre located at (awew,0). We ac- 
count for the position of periastron (defined in figure 1) by 
rotating the ellipse counter-clockwise by an angle vow about 
the z'-axis. This gives us the S[ frame: 



Xi — x cos zuw — y smzvw 
y[ — x sin vjw + y cos zuw 



(A3) 
(A4) 



To match the figure, we require a translation to place 
the planet-moon barycentre at the origin. After applying 
this translation, we have found the desired S' 2 frame. 



X2 = X\ — aw&w cos zuw 
2/2 = Vi ~ a w ew sinzuw 



(A5) 
(A6) 
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Note that ew, the eccentricity of the planet's orbit 
around the planet-moon barycentre, must be equal to es- 
However, zuw will have a phase difference of it relative to 
vjs- In the defined ellipse, the distance between the barycen- 
tre and the centre of the planet, rw as a function of true 
anomaly, fw, is given by: 



rw(fw) = 



Ow(l 



e 2 w ) 



(A7) 



1 + ew cos fw 

With this equation, we may define the x' and y' position 
of the centre of the planet as a function of fw ■ 

x' = a w ew + r w (fw) ■ cos fw (A8) 
V = rw(fw) ■ sin f w (A9) 

Having defined x'(fw) and y'(fw), we have also ac- 
quired x' 2 (fw), using equation (A5). Putting x' 2 (fw) into 
equation (A2) and integrating between the limits and 2n, 
as dictated by (Al), we find: 



s aw ■ Cries, ^s) 

Ottv = —j= - 

V2 ■ v B j_(e P , zu P ) 

where we define (ries^s) by: 



(A10) 



(l-e|) 1/4 



es 



\Je 2 s + cos(2o7 5 )(2(l - e|) 3 / 2 - 2 + 3e 



(All) 

Note, that we have replaced ew with es, since they are 
equivalent, and zuw with zus since the n phase difference 
between them does not affect this expression. 

aw may be re-written as: 

Ms 

a w = a s — (A12) 

Mprv 

And as may be expressed in terms of the Hill radius, 



as = xdn 



xap {-^m — ) 



(A13) 



In this expression, ap may also be re- written using Ke- 
pler's Third Law: 



ap 



-( 



G(Mp RV + M,)Pp \ i/3 
4tt 2 / 



(A14) 



Finally giving: 



aw 



X ■ 



Ms (G{M PRV +M*){M PRV - M s )Pp\ 1 / s 



M F 



12tt 2 M, 



(A15) 

Let us now turn our attention to the expression for 
VB±(ep,vjp). If we assume the time it takes for the planet 
to cross a distance aw is small compared to the period of 
the orbit, then we may assume vb± does not vary. However, 
it is only the perpendicular component of vp , which we label 
as vp±, which dictates the TTV effect. In the case of a cir- 
cular orbit, it is easy to see that vp± = vb, but for eccentric 
obits, this is not the case. 

Consider the orbit of the planet-moon barycentre 



around the host star. To find the perpendicular component 
of the barycentre's orbital velocity, we must find the gradi- 
ent of the tangent to an ellipse which has been rotated for 
the position of periastron, zup, and then translated so that 
the star is at the origin. 

We now consider an x-y plot in a frame S with an ellipse 
centered at the origin and the star located at coordinates 
(apep, 0). We then apply a counter-clockwise rotation about 
the 5-axis by an angle zup into the Si frame. We then apply 
the translation to reach the & frame. The equations are 
identical to (A3) & (A4) and (A5) & (A6) except we replace 
the W subscript with P. These equations may be re-written 
to make x & y the subject as follows: 



x = —x~2 cos zup — y~2 sin vjp — apep (A16) 
y = x~2 sinrop — y~2 coszup (A17) 

Note that we have assumed the orbital eccentricity of 
the planet-moon barycentre around the host star is equiva- 
lent to the orbital eccentricity of the planet around the host 
star. Equations (A16) & (A17) must satisfy the standard 
equation of an ellipse: 



+ ■ 



(A18) 



We substitute equations (A16) & (A17) into equation 
(A18) and then differentiate implicitly with respect to x~2- 
We then rearrange to make &y~2/&X2 the subject. Finally, we 
substitute for £2 and y 2 using equations (A5) & (A6) to get: 



dj/2 

dX2 



y sin-cup — (1 



e P )£ coszup 



y coszup + (1 — e 2 p )xsinzup 



(A19) 



x(fp) and y(fp) are known through the equations (A7), 
(A8) and (A9) except the subscript W should be replaced 
by P in these three expressions. 

The angle the gradient of the tangent makes to the +x 
axis is given by: 



0(fp) = arctan 



[ dy~2 ] 
I. deb J 



(A20) 



And finally we may express the perpendicular compo- 
nent of the barycentric velocity as: 



v B ± = v B cos[0(/p)] 



(A21) 



fp is the true anomaly of the planet-moon barycentre 
during transit and is given by fp = it /2 — zup. Performing 
the necessary substitutions we find: 



«»*(/„)= cos [arctan (^^)\ (A22) 

If we take the limit of the expression (A24) in the case of 
ep — > 0, we get the expected result of simply vb- vb is also a 
function of ep and zup and there is an obvious desideratum 
to create a single parameter which absorbs all the effects of 
planetary eccentricity. Thus we write % as: 
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Vb 



G(M PRV + M,) 2(1 + e P sinrop) _ 1 



(A23) 



And we may now create a single parameter to absorb 
the effects of planetary eccentricity, T: 



~- r / — ep cos vjp \\ 

T = cos arctan 

L V 1 + ep smvop / J 



epcostup y\ 1 2(1 + ep sin tup) 



- i 



(A24) 



Thus: 



tv ^ G(Mprv + M,) 
v B ± = T(e P ,wp) ■ \ — (A25) 



fflp 



The final expression may be written as: 



5 TTV = aw ^ — . ^s,™s) (A26) 

V2^G(M PRV +M*) T(ep,uj P ) 

where aw is given by equation (A15), ap is given by 
equation (A14), £ T (es,vjs) is given by equation (All), 
T(ep,wp) is given by equation (A27). An alternative ex- 
pression of (A26) is given by: 



1 a y, 2 a s M s Mp RV C, T ( 

OTTV — — 1= ■ — = ■ ttt 7 (AZl ) 

V2 \J G(M* + Mprv) T(ep,rop) 

It is therefore clear that our derived equation is equiva- 
lent to that of previous authors in the case of circular orbits. 
It is also clear that Sttv oc Msas- 

For completion and scaling purposes, we also derive the 
TTV amplitude in terms of just masses and periods by using 
Kepler's Third Laws. 



In comparison, the orbital period of the planet around 
the host star is given by: 



2„3 
P 



,G(M,+M PRV ] 
Therefore we may write: 

Ps_ „ 
Pp ~ 



1/2 ^ /47r 2 apy/2 
~ V GM, 



(B4) 



(B5) 



APPENDIX C: DERIVATION OF THE TDV 
AMPLITUDE 

CI Derivation of the circular form 

We will follow the same definition for TDV as for TTV, 
where TDV is equal to the observed - calculated (O-C) tran- 
sit duration. The calculated (or expected) transit duration, 
r, is given by: 



X 

vb± 



(CI) 



Where X is the distance the planet has to cross in order 
to complete the transit and vb is the velocity of the planet 
during transit, as given by Kepler's Third Law. 

In the case of an additional moon orbiting the planet, 
the velocity of the planet now has two components, the ve- 
locity of the planet-moon barycentre (vb) and the wobble 
velocity due to the perturbation of the moon. We remind 
the reader that we assume coplanar orbits with i = 90° . 



r(fw) 



X 



Vb± + vw±(fw) 
Therefore, TDV is defined as: 



(C2) 



ottv = — • - - 7= • Z T (M,,M PR v, Ms) ■ -=r, r 

2n S 1 /- 1 ^ T(ep,wp) 

(A28) 

where Z T (M t , M PRV , M s ) is the TTV mass function 
given by: 



/ Mf(AW-M s ) 2 y /6 ,, ( Ml y/e 
T ~\ M% RV Ml ) -\M2M PRV ) lA29j 

Where the approximation is based on Mp R v 3> Ms- 



TDV(fw) =f- r(f w ) = ( 

V vb± + 



Vbj_ 



■ v w ±(fw) 
If v B ± 3> v w ± , then we may write: 



TDV(/ W 



vw±(fw) - 

• T 

V B ± 



l)-f (C3) 



(C4) 



If we assume a circular orbit, then the velocity of the 
planet around the barycentre of the planet-moon system, 
the wobble velocity, is given by: 



APPENDIX B: THE P S /P P RATIO 
From Kepler's Third Law, we may write: 



2nav 



v w = 



Ps 



2TvM s a s 1 
Mprv ' ~P~s 



Using equation (Bl), we may write: 



(C5) 



\GM PRV ) 

where as is some fraction, x, of the Hill radius: 



(Bl) 



as = X-a P - 1 



M PRV - M s \ 



3M* 



(M PRV \ i/3 

ap \^M7) (B2) 



X-a P - 



Putting (Bl) and (B2) together: 

'47rVa 3 P \i/2 



b \ 3GM» J 



(B3) 



vw 



asMp RV 



(C6) 



However, we want the perpendicular component of this 
velocity, v w ± as a function of true anomaly, fw- For a cir- 
cular orbit, it is trivial to show that the variation will be 
sinusoidal and hence the rms amplitude of v w ± is given by: 



\Vw±\ 



1 / GM's y/ 2 
s/2 \a s Mp RV ) 



(C7) 
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Let us now consider the velocity of the planet-moon 
barycentre around the host star. For a circular orbit, the 
velocity of the barycentre is given by: 



Vb = V B A 



2-wap 



(C8) 



Pp V P P 

Once again, due to the circular orbit, the perpendicu- 
lar component of the velocity is the same as the absolute 
value during transit. The ratio of v w ±/vb± multiplied by 
the transit duration will be the TDV signal. 



TDV(f w ) = - Vw±{fw) -f = -R.f (C9) 
vb± 

In the case of circular orbits, R is therefore given by: 



the perpendicular component of the planet's wobble velocity 
to the perpendicular component of the velocity of the planet- 
moon barycentre around the host star, multiplied by the 
duration of the transit. 



TBV(f w ) = vw,±(M 

Vb± 



(C17) 



In order to find the rms amplitude of this signal, wc 
need to derive vw,±(fw), which is not the same as vw(fw)- 
To do this, we need to find the gradient of the tangent of the 
planet's position along our transformed ellipse from figure 
A. To find the gradient of the tangent to the ellipse at any 
true anomaly, fw, we need to implicitly differentiate the 
equation for the ellipse in the S' 2 frame. Equations (A5) and 
(A6) may be rewritten making x and y the subject: 



R _ 1 / GM 2 S y/2/ 27rGM. \-V3 



^asMpavJ < - ' (C10) 

From equation (C5), we may instantly infer that the 
TDV effect is oc M s ag 1/2 , in contrast to the TTV effect 
which is oc Msas. In the case of M„ 3> Mppv and Mppv >• 
Ms, as is given by: 



it i 
x — —x 2 cos ujw — J/2 sin ujw — awew 

y — x 2 sin ww — y'i cos zuw 



(CIS) 
(C19) 



These two equations must satisfy the standard equation 
of the ellipse: 



(GMppvPp^/ 3 
And feeding this into the equation for R we get: 



R = J_ J_ ( 3M| 



1/6 



^2 VX \MtM PRV < 



(Cll) 



(C12) 



Therefore the TDV amplitude for a planet, on a circular 
orbit, with a single moon, which is also on a circular orbit, 
is given by: 



t 3 1/6 

Stdv ^ — • ■ Z D (M,,M PRV ,Ms) 
VX V2 

where the TDV mass function is defined by: 



Z D {M,,Mp RV ,Ms) ~ ( 



M% 



1/6 



\M3M PRV , 



(C13) 



(C14) 



+ 



(C20) 



We substitute (C18) and (C19) into (C20) and then 
differentiate implicitly with respect to x' 2 . We arrange the 
equation to make dy 2 /dx 2 the subject. Finally, we replace 
the x' 2 and y' 2 terms using equations (A5) and (A6) and find: 



dj/2 y' smzow ~ (1 — e^x' coszaw 
dx' 2 y' coszuw + (1 — £w) x ' sinwiy 



(C21) 



x' (fw) and y'(fw) are known and so we have found the 
gradient of our tangent. The angle of the tangent is given 
by: 



9'(h 



arctan 



[Ml 



(C22) 



And finally we may express vw,±(fw) as: 



C2 Derivation of the eccentric form 

Here, we derive the TDV rms amplitude in the case of eccen- 
tric orbits. We use the definition of TDV given by equation 
(C4). The orbital velocity of the moon around the barycen- 
tre of the planet-moon system for a circular orbit is given 
by: 



vw,±(fw) = vw(fw) ■ cos[6'(fw)] 



(C23) 



We now take the rms of this waveform in the way de- 
tailed in appendix A and find: 



Stdv 



t 

vb± 



IMV 
2aw 



(C24) 



v w = 



GMg 
asMpRv 



GM% 



(C15) 



In the case of an eccentric orbit, this velocity becomes 
a function of true anomaly, fw ■ 



vw(fw 



1/2 

Vw 



-4 

aw / 



1/2 



(C16) 



<r w (fw) a w , 
The TDV signal has a waveform governed by the ratio of 



where 



Co(es, ru s ) 



1 + e| — e| cos(2tj7s) 



i-4 



(C25) 



Note that we have replaced ew by es since they are 
equivalent and similarly for cos2zow and cos 2ms. Taking 
the limit of this expression as es — > gives the expected 
result of Ty / nw/vB\/2aw- v B ± has already been derived in 
appendix A and equation (A28). 
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~, s /G(M PRV + M t ) 

v B ± = i{e P ,wp) ■ \ — (C26) 

V ap 

Using this and equation (C9) we may re-write Stdv as: 



a P Ml t ( D (e s ,KJs) 



as y Mprv(Mprv + Af*) ^2 T(ep,ro P ) 

(C27) 

This once again demonstrates the eflect's proportional- 
ity of oc Msa^ 1 ' 2 . Using equation (A13), the fraction ap/as 
may be substituted for and we find: 

Stdv = ~^-Z D (M ry Mp RV ,M s y C ^ eS ' WS } (C28) 
VX V2 T(ep,zu P ) 

where 

M„Mf 



Z d {M*,Mprv,M s ) = (- 



,Ml RV {M PRV + M*y{M PRV - M s ), 

(C29) 

This general expression can be shown to be equiva- 
lent to our approximate form with the same approximations 
made. 



APPENDIX D: RATIO OF TDV TO TTV 
AMPLITUDE 

In this section we derive the ratio of the TDV amplitude 
to the TTV amplitude. This will allow us to analytically 
quickly see which signal is stronger for any new system we 
come across and solve for the exomoon mass exactly. We 
have: 



PP X Z T (M„Mp RV ,Ms) Ct(6S '^ 



2tt 3V3V2 V ' ' T(e P , rop ) 

(Dl) 



Stdv — —7= ■ —7=- ■ Zd (M, , M PRV , M s ) ■ ijffr ' Ws } (D2) 
y/X V2 T(ep,ra7p) 

We make the assumptions that M* 3> Mp R v and 
Mp R v ^> Ms- In this case, the TTV mass function, Zt, 
and the TDV mass function, Zd, are equal. We also make 
the assumption that moon is on a circular orbit and so 
Ct = Co = 1- Furthermore, the planet's eccentricity pa- 
rameter, T will cancel out and so we need not make the 
assumption ep = 0. This gives us the ratio of TDV to TTV, 
ij to be: 

_ Stdv 27tt \/3 _ 2itf . . 

v - w " Tp~ ' W ~ ~K ( ' 

The masses have completely cancelled out and we are 
able to solve the equation for Ps in a completely rigourous 
way. 
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by the author. 



